1 



Lecture Notes on Fluid Dynamics 

(1.63J/2.21J) 
by Chiang C. Mei, 2002 



7.3 The Shallow- Water Approximation 

The continuity equation reads 

du ^ dv ^ dw Q 

dx dy dz 



u ^ 



If the horizontal length scale of motion L is much greater than the sea depth D, i.e., 
shallow water, 

^ < 1 (7.3.1) 



we have 



Momentum conservation requires 



L 

W^^^U (7.3.2) 



du ^. du ^ du ^ du ^ 1 dpa 

T L — L pL 

dv ^ ( dv ^ dv ^ 9v\ ^ ^ 1 dpd 
dt \ dx dy dz J P dy ^ 

dw ^ , dw ^ dw ^ dw . 1 dpa 

^dt^ ^_dx^ dy ^_dz, p dz 

T L — L 

where Pd stands for the dynamic pressure 

p = -pgz + pd 

Examine first the horizontal momentum balance. In order that pressure gradient forces 
the rotating flow, the pressure gradient has to balance either the Coriolis force (rotation 
important), assuming that the Rossby number is no greater than 0{1). The scale of pressure 
is 

P = ^, or ,UfL 



Prom the vertical momentum we find, if P = pUL/T, 



dw DU 7-)2 

dt LT ^ 



Idd pUL/T 12 
p dz pD 
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and if P = pUfL, 

dw DU jj2 ^ 



1^ ML L2 Tf 

p oz pD ■' 

Since the time scale of interest is of a day or so, /T = 0(1). We conclude that the vertical 
pressure gradient is practially zero with error of order D'^ jl? <^ 1, implying that 

Vd ^ pgr] (7.3.3) 

or the total pressure is hydrostatic 

^^total - ^^(^ - ^) (7-3-4) 

Thus, 

dvd df) dpd Of] 

We now show by a formal perturbation scheme for horizontal bottom, that u, v are 
independent of z to the leading order or approximation, 

du dv Q (73 g^ 
dx dy dz 

— + u— + v— + w—-fv = -— (73 7) 
dt ^ dx ^ dy ^ dz ^ ^ dx 

dv dv dv dv , dn ,^ . 

w^O (7.3.9) 



On the bottom 



On the free surface 



df) dri dri ^701^ 



Assuming the following expansions: 



z + h {z + hf /^o-.-.N 

u^uo + -J— ui+ ^ U2 + --- (7.3.11) 

v = vo + -J— VI+ '- V2 + --- (7.3.12) 

z + h , {z + hy 

w = wi H ^J2~ ^2 H (7.3. 13) 
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where Un,Vn, Wn are independent of z. Prom the x momentum equation we get from (7.3.7) 

duQ , duQ , duQ . , dr] 

— + Uq—- + Vq— Jvq + g— 

at ox ay ox 

(z + h) ( dui duo dui duo dui wiUi \ 
+ ui— \-uo^ ^vi— \-vo— \ JVi] 



L \ dt dx dx dy dy 

+ (^^Ji...) + ... = (7.3.14) 

with a similar equation for the y momentum. Separating the zeroth power of {z + h)/L we 
get 

duo duo duo dr] ^7 o i r x 

dvo dvo . -.dvo dr] 

^ + ^.o^ + .0)^ + fuo = -9g^ (7.3.16) 

Thus u, V are depth-independent to the leading order. Also from continuity, 

dun dvn Wi „x 

and from the free surface condition 



dn dri ri + h /^, o 



Hence after eliminating wi, 



dv ^ d[{7] + h)uo\ ^ d[{7] + h)vo\ ^ ^ 
dt dx dy 

In summary, for shallow seas, 

du du du , dn ^ 

dt dx dy dx 

dv dv dv „ dn ,^ . 

dt dx dy dy 

to leading order. 

If Rossby number is small, then the convective inertia is negligible; the momentum equa- 
tions reduce to : 

where are linear. 

Take for estimates, U = 0.1,lm/s,n = 2.31 x 10-^L = 100, 1000A;m = 10^ lO^m, 
Rossby number U/2QL = 0.043. We leave it as an exercise to work out the equations for 
h{x, y) with a small slope. 
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7.3.1 Geostrophic motion 

For steady flow at small Rossby number, 

H — rj + h c:i h 

the momentum equations reduce to 

Thus, Coriolis force and pressure gradient are in balance 

g dr) g dr] 

J dy J dx 

implying 

«^ + -^ = (7.3.24) 

or 

q- V?7 = 0. 

Physically along a streamlines, the free surface height remains constant. Hence, the sur- 
face contours are parallel to the streamlines and to isobars. This state is called geostrophic. 

7.3.2 Remark: Depth-integrated mass conservation 

The depth integrated mass conservation (7.3.25) is a general an exact result which holds for 
variable depth as well. On the seabed z — —h{x, y), vanishing of the normal veocity requires 

dh dh 
w = -u- V—. 

ox oy 

On the free surface z — r]{x, y, t), the kinematic boundary condition reads, 

dr] dr] dr] 
w = — + M— + u— z = rj 
at ox oy 

Integrating the continuity equation and using Leibniz's rule 

n ( du dv dw\ , 

r n») dp, dp, 

^ ^''^--^Yxl-n'''^'^dylV'^' 

dr] dr] dh dh 
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Define the total depth tohe H — r] + h 

^ ( dr] dri\ ( dh dh\ dHu diHv) 
u—\w — u— v—\ —\w + u— — h f - 



dx j \ J -h ^-^ '^y 

Using the boundary conditions we get 

dH dHu dHv _ ^ (7 3 25) 

dt dx dy 



